arXiv:1502.05889vl [nucl-th] 20 Feb 2015 


Gamow-Teller properties of the double /3-decay partners ^^®Cd(Sn) and ^^°Nd(Sm) 


D. Navas-Nicolas and P. Sarriguren 

Instituto de Estructura de la Materia, lEM-CSIC, Serrano 123, E-28006 Madrid, Spain 

(Dated: February 23, 2015) 

The two Gamow-Teller (GT) branches connecting the double-beta decay partners (^^®Gd, ^^®Sn) 
and ('^^^Nd, with the intermediate nuclei and ^®°Pm are studied within a microscopic 

approach based on a deformed proton-neutron quasiparticle random-phase approximation built on 
a Skyrme selfconsistent mean field with pairing correlations and spin-isospin residual forces. The re¬ 
sults are compared with the experimental GT strength distributions extracted from charge-exchange 
reactions. Combining the two branches, the nuclear matrix elements for the two-neutrino double¬ 
beta decay are evaluated and compared to experimental values derived from the measured half-lives. 
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I. INTRODUCTION 

Neutrinoless double-beta decay {QvjSjS) is nowadays a 
topical issue that has attracted a lot of interest from both 
theoretical and experimental sides [H. This decay mode 
is still unobserved and would violate lepton number con¬ 
servation. Its existence would demonstrate that the neu¬ 
trino is a massive Majorana particle and would provide a 
measurement of the absolute mass scale of the neutrino. 
Obviously, the nuclear matrix element (NME) involved in 
such a process must be determined accurately to extract 
a reliable estimate of the neutrino mass. On the other 
hand, the two-neutrino double-beta decay mode {2uj3(3) 
is perfectly allowed by the Standard Model. It is a rare 
second-order weak-interaction process that has been ob¬ 
served experimentally in several nuclei (see Ref. 0 for 
a review). Thus, to test the reliability of the nuclear 
structure calculations involved in the Quj3(3 process, one 
checks first the ability of the nuclear models to reproduce 
the experimental information available. 

Certainly, there are differences between the NMEs in¬ 
volved in both processes, but there are also clear similar¬ 
ities. Among the differences, one should notice that the 
2^/3/3-decay NMEs are dominated by Gamow-Teller (GT) 
transitions connecting the initial and final 0+ states with 
all J'^ = 1+ states in the intermediate nucleus. On the 
other hand, in the Oiz/3/3 process the intermediate states 
run over all J'^ values. Another difference with conse¬ 
quences in the evaluation of the NMEs is that in 0^/3/3 
decay the average virtual neutrino momentum is much 
larger than the typical scale of nuclear excitations and 
then closure approximation can be safely used. In the 
case of 2vf3P decay the neutrino momenta are compara¬ 
ble with the nuclear excitation energies, preventing the 
use of closure. Among the similarities, the two processes 
connect the same initial and final nuclear ground states 
to be described within a given nuclear model and they 
share the calculation for the intermediate J'^ = states. 
Therefore, reproducing the 2vPP NMEs is a requirement 
for any nuclear structure model aiming to describe the 
neutrinoless mode. One can go even further and compare 
the calculations, not only with the 2vfij3 NMEs extracted 
from the measured half-lives, which is nothing but a num¬ 


ber, but also with the GT~ (GT'*') strength distribution 
of the single branch connecting the initial (final) ground 
state with all the J'^ = I’*' states in the intermediate nu¬ 
cleus. This comparison provides a more detailed informa¬ 
tion on which the calculation can be tested. GT strength 
distributions have been measured in recent years from 
high resolution charge-exchange reactions (GER), within 
a large experimental program aimed to explore the GT 
properties at low excitation energies of double /3-decay 
partners M- 

The purpose of this work is to investigate the possibil¬ 
ity to describe the rich information available at present 
on the GT nuclear response within a formalism based on 
a deformed proton-neutron quasiparticle random-phase 
approximation (QRPA). This information includes global 
properties such as the location and strength of the GT 
resonance, more detailed description in the low-lying ex¬ 
citations, and last but not least, the 2:z/3,5-decay NMEs. 
This study was started for the double /3-decay partners 
with A = 76,128,130 d, d|. Here, we shall focus on 
the ii®Gd ^ ii®Sn and ^'^“Nd ^ i5°Sm decays, moti¬ 
vated by the recent experiments on GER performed for 
^^®Gd(p,n)^^®In and ^^®Sn(n,p)^^®In Q, as well as for 
i50Nd(3He,t)i'50pm and i50Sm(t,3He)i5°Pm 0. It is also 
worth mentioning the recent measurement of the electron 
capture in ^^®In [l^ that can be used as a benchmark for 
nuclear structure calculations. In addition, the double /3 
decay of ^®°Nd has received increasing attention in the 
last years and it is currently considered as one of the best 
candidates to search for the OuP/d decay in the planned 
experiments SNO-I-, SuperNEMO, and DGBA. The rea¬ 
son for the interest in ^^°Nd is that it has a large phase- 
space factor and therefore a relatively short half-life. It 
also has a large Qpp energy that will reduce the back¬ 
ground contamination. However, both ^^°Nd and ^®°Sm 
are deformed nuclei that require a deformed formalism 
to deal with them properly. 

The paper is organized as follows: In Sec. m we 
present a brief summary of the theoretical approach used 
to describe the GT strength distributions, as well as the 
2^/3/3-decay basic expressions. SectionHIIlcontains the re¬ 
sults obtained from our approach, which are compared to 
experimental data and other available calculations. Sec- 
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tion m contains a summary and the main conclusions. 


II. THEORETICAL APPROACH 

The description of the deformed QRPA approach used 
in this work is given elsewhere [111411^ . Here we give 
only a summary of the method. We start from a self- 
consistent deformed Hartree-Fock (HF) calculation with 
effective nucleon-nucleon Skyrme interactions, assuming 
axial deformation and time reversal symmetry [l^ . Most 
of the results in this work are performed with the Skyrme 
force SLy4 [l^ , which is one of the most widely used and 
successful interactions. In addition, we also present re¬ 
sults obtained with other Skyrme forces to investigate 
their suitability for the description of the spin-isospin 
properties of nuclei. In particular, besides SLy4, we shall 
show results from the simpler force Sk3 [l6l| and from 
the force SGII [13 that was fitted taking into account 
nuclear spin-isospin properties. 

In our approach, the single-particle wave functions are 
expanded in terms of the eigenstates of an axially sym¬ 
metric harmonic oscillator in cylindrical coordinates us¬ 
ing twelve major shells. Pairing correlations between 
like nucleons are included in BCS approximation taking 
fixed pairing gap parameters for protons and neutrons, 
which are phenomenologically determined. Besides the 
selfconsistent HF solution, we also explore the potential 
energy curves, that is, the HF energy as a function of 
the quadrupole deformation /3, which are obtained from 
constrained HF-I-BCS calculations. These results can be 
seen in Fig. [T] for ^^®Cd and ^^®Sn, and in Fig. [5] for 
^®°Nd and ^^°Sm for the three Skyrme forces considered. 
We get very soft profiles for ^^®Cd and ^^®Sn, whereas for 
^®°Nd and ^^°Sm we obtain two energy minima, oblate 
and prolate, but with clear prolate ground states in both 
cases. We obtain similar results with the three Skyrme 
forces with the only exception of the force Sk3 in ^^®Cd, 
which seems to miss a spherical solution. 
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FIG. 1: (Color online) Energy-deformation curves for 
(a) and ^^®Sn (b) obtained from HF-I-BCS calculations with 
various Skyrme forces. 
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FIG. 2: (Color online) Same as in Fig. [T] but for ^®°Nd (a) 
and i®“Sm (b). 
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FIG. 3: (Color online) GT strength distributions B{GT~) in 
'^^^Cd as a function of the excitation energy in the daugh¬ 
ter nucleus (a) and accumulated strengths (b). SLy4-QRPA 
calculations for various deformations are compared to data 
(Sasano 2012: [^) from (p, n) reactions. 


After the HF-I-BCS calculation, the QRPA equations 
are solved on the deformed ground-state basis to get the 
GT strength distributions and to compute the 2vf}j3- 
decay NME. To describe GT excitations in QRPA we add 
to the quasiparticle mean field a separable spin-isospin 
residual interaction in the particle-hole {ph) and particle- 
particle {pp) channels. The ph part is responsible for the 
position and structure of the GT resonance and its cou- 
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FIG. 4: (Color online) GT strength distributions B(GT’'') in 
^'^®Sn as a function of the excitation energy in the daugh¬ 
ter nucleus (a) and accumulated strengths (b). SLy4-QRPA 
calculations for various deformations are compared to data 
(Sasano 2012: [^) from (n,p) reactions, as well as from 

(d,^He) (Rakers 2005: i). 


pling constant Xph usually taken to reproduce them 
[H The pp part consists of a proton-neutron 

pairing force and it is also introduced as a separable force 
[l^ . The coupling constant is usually fitted to 
the half-lives phenomenology [2^. In this work we have 
chosen for the ph and pp coupling constants the law de¬ 
pendence on the mass number A given by C in Ref. 

However, the different mean fields used, Nilsson 
20j versus Skyrme-HF here, make the constants some¬ 
what different. Specifically we use Xph — 3.0/H° ’^MeV 
and MeV. Although slightly different val¬ 

ues could be used for the different Skyrme forces, we use 
in this work the same values for them, so that changes 
among calculations from various forces are caused by the 
underlying quasiparticle structure provided by the effec¬ 
tive Skyrme force. At this point it is worth mention- 
ing the studies on ^®°Nd and ^®°Sm carried out in Refs. 
21, where deformed QRPA calculations were per¬ 
formed using realistic nucleon-nucleon residual interac¬ 
tions based on the Brueckner G matrix for the CD-Bonn 
force on top of a phenomenological deformed Woods- 
Saxon potential. Comparison with the results obtained 
from schematic separable forces [2l| shows that the latter 
reproduce fairly well the main characteristic of the GT 
strength distributions profiles and 2i//3,0-decay NMEs. 
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FIG. 5: (Golor online) Same as in Fig. [3l but for B(GT ) in 
i®°Nd. Data (Guess 2011: 01) are from (®He,t) CERs. 


The technical details to solve the QRPA equations have 
been described in Refs. For each value of the 

excitation energy, the GT transition amplitudes in the 
intrinsic frame connecting the ground state |0) to one- 
phonon states in the daughter nucleus \ujk), are found to 
be 

{uJK\<JKt^\0)=TM^^ , ( 1 ) 

where 

TTt' 

TTt' 

in terms of the occupation amplitudes for neutrons and 
protons (ut tt = 1 ~ the forward and back¬ 

ward amplitudes of the QRPA phonon operator, 
and , respectively. Once the intrinsic amplitudes in 
Eq. (P) are calculated, the GT strength B(GT) in the lab¬ 
oratory frame for a transition IiKi{0'^0) —>■ IfKf{l'^K) 
can be obtained as 


B<^(GT=^) = V [cTot^l O)^ ()(a;jc^o -uj) 

(jJK 

- 1-2 {ujk=i I O )^ 5(ijJk=i — w ) 


■(4) 
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FIG. 6: (Color online) Same as in Fig. H] but for B(GT''') in 
Data (Guess 2011: Q) are from (t,®He) CERs. 


To obtain this expression we have used the initial and 
final states in the laboratory frame expressed in terms 
of the intrinsic states using the Bohr and Mottelson fac¬ 
torization. Finally, a quenching factor q = gAlgA,hai:e = 
0.79 is included in the calculations to take into account 
in an effective way all the correlations that are not prop¬ 
erly considered in the present approach. We shall discuss 
later in Subsec. mm various attempts aiming to fit si¬ 
multaneously single-,5 and double-/3 decay observables by 
adjusting the quenching factor. 

Concerning the 2i//3/3-decay NMEs, the basic expres¬ 
sions for this process, within the deformed QRPA formal¬ 
ism used in this work, can be found in Refs. [T^. [^ [^. 
Deformation effects on the 2z//3/3 NMEs have been also 
studied within the Projected Hartree-Fock-Bogoliubov 
model [ 2 ^. Attempts to describe deformation effects on 
the Ovpp decay within QRPA models can also be found 
in Refs. [ 2 ^ [ 23 . The half-life of the 2i//3/3 decay can be 
written as 



FIG. 7: (Color online) GT strength distributions in (a), 

^^®Sn (b), ^®°Nd (c), and ^®'^Sm (d). The results correspond 
to the prolate shapes with the Skyrme forces SLy4, SGII, and 
Sk3. Data are as in the previous figures. 


rTi2,iy 0 j3 
^ 1/2 



= (gA)^ 





in terms of the phase-space integral and the nu¬ 
clear matrix element that contains all the infor¬ 

mation of the nuclear structure involved in the process, 


mi,mf 

x{0f\a-Kt~\uJK,mf) (Wiy,mi |0i) . (6) 


In this equation \u}x,mi){\^K,mf)) are the QRPA inter¬ 
mediate I’*' states reached from the initial (final) nucleus. 
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FIG. 8: (Color online) Same as in Fig. [S] but in the low 
energy range. Also shown is the GT str eng th extracted from 
electron capture of ^^®ln (Wrede 2013: [Q)- 


FIG. 9: (Color online) Same as in Fig. [4] but in the low 
energy range. Also shown (open circle) is the GT strength 
extracted from f5~ decay of ^^®In to which the data from 
are normalized. 


The indices rrii, mf label the 1+ states of the interme¬ 
diate nucleus. The GT matrix elements are those in 
Eq. ([T]). The overlaps are needed to take into account 
the non-orthogonality of the intermediate states reached 
from different initial | 0 i) and final | 0 f) ground states. 
Their expressions can be found in Ref. (23|. The energy 
denominator D involves the energy of the emitted lep¬ 
tons, which is given on average by \Qpp + We, as well 
as the excitation energies of the intermediate nucleus. In 
terms of the QRPA excitation energies the denominator 
can be written as 

Z?i = i(a;r+<-), (7) 

where ) is the QRPA excitation energy relative 

to the initial (hnal) nucleus. It turns out that the NMEs 
are quite sensitive to the values of the denominator, espe¬ 
cially for low-lying states when the denominator reaches 
smaller values. Thus, it is a common practice to use some 
experimental normalization of this denominator trying to 
improve the accuracy of the NMEs. In this work we shall 
also consider the denominator Z? 2 , which is corrected 
with the experimental energy of the first 1 '*" state 
in the intermediate nucleus relative to the mean ground- 
state energy of the initial and hnal nuclei, in such a way 
that the hrst calculated state appears at the exact 
experimental energy. 


D2 
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+ LO 


rrii 
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+ + . 
-*-1 


( 8 ) 


The running 2://3/3 sums calculated later, will be shown 
for the two choices of the denominator Di and 02- In 
the case of the A = 116 isobars, the ground state in 
the intermediate nucleus ^^®In is a I’*' state. Then, the 
energy is given by 


^1+ — 2 T Qf}- )exp ) ( 9 ) 

written in terms of the experimental energies Qec = 
0.463 MeV and = 3.276 MeV of the decays of 
^^®In into ^^®Cd and ^^®Sn, respectively. In the case of 
A = 150, although the ground state in ^^°Pm is not a 
state, the hrst I'*' state excited in CERs appears at 
a very low excitation energy, E=0.11 MeV Q. Hence, 
we also determine from the experimental values of 
Qec = 0.083 MeV and Qa- = 3.456 MeV for the decays 
of ^®°Pm into ^®°Nd and ^^°Sm, respectively. 

The various measurements reported for the 2v(3j3- 
decay half-lives have been analyzed in Ref. [1], 

where recommended values of r^^y^ 2 ^^(^^®Cd)=( 2 . 8 ± 0 . 2 ) x 
10^9 yr and Q^°m)={8.2 ± 0.9) x 10^® yr were 
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FIG. 10: (Color online) Same as in Fig. [5l but in the low 
energy range. 


adopted. Using the phase-space factors from the re¬ 
cent evaluation (28j that involves exact Dirac wave func¬ 
tions with finite nuclear size and electron screening, 
Q 2 yfip (ii6cd)=2.76410-1® yr-i, and (i®0Nd)= 

3.64310-11’ yr-i, one gets the experimental nuclear 
matrix elements (ii®Cd)=0.138 MeV-i (0.223 

MeV-i) when the bare gA,hare = 1.273 (quenched fac¬ 
tor gA = 1) is used, and (i®°Nd)=0.070 MeV-i 

(0.114 MeV-i) when ^^.bare = 1.273 {gA = 1) is used. 


III. RESULTS 

A. Gamow-Teller strength distributions 

Next we discuss the results obtained for the energy dis¬ 
tributions of the GT strength. In the upper panels (a) 
of Figs. [3] and S] we can see the GT strength distribu¬ 
tions for ii®Cd (i3(GT-)) and n^Sn (i3(GT''')), respec¬ 
tively. The lower panels (b) show the respective accumu¬ 
lated GT strengths. We compare our QRPA results from 
SLy4 obtained for prolate, oblate, and spherical shapes 
in the initial (H^Cd) and final (H^Sn) nuclei with the 
experimental strengths extracted from ii®Sn((i,^He)ii®In 
CERs 0 and from ii®Cd(p, n)ii®In and ii®Sn(n,p)ii®In 
CERs Q . One should notice that the measured strength 
extracted from the cross sections contains two types of 



E [MeV] 

FIG. 11: (Golor online) Same as in Fig. [6] but in the low 
energy range. 


contributions, namely GT {at^ operator) and isovector 
spin monopole (IVSM) {r^at^ operator). Both are as¬ 
sociated to AL = 0 and AS' = 1 transitions and there¬ 
fore, they have similar angular distributions that can¬ 
not be disentangled in the experiment. Thus, the mea¬ 
sured strength corresponds actually to B(GT-flVSM). 
Different theoretical calculations evaluating the contri¬ 
butions from both GT and IVSM modes are available 
in the literature [1, [2§-[3lj . In Ref. a self-consistent 
Skyrme Hartree-Fock plus Tamm-Dancoff approximation 
was used to separate the GT from the IVSM strength. 
Well separated centroids of both resonances were found, 
with the IVSM resonance at much higher energy, as it 
corresponds to a 2hui mode. Subsequent calculations in¬ 
cluding A = 116 isobars were performed in Refs. 0,1111 
with different degrees of complexity, from a schematic 
single-particle model space up to realistic calculations 
within a QRPA approach in a Woods-Saxon basis and 
two-body interactions constructed with the Bonn-A po¬ 
tential. In Ref. @ the contribution from the IVSM 
component was also evaluated by employing a micro¬ 
scopic method based on QRPA with residual interac¬ 
tions obtained from the Brueckner G-matrix for CD- 
Bonn nucleon-nucleon force. In general, the conclusions 
of those theoretical calculations converge in an overall 
picture. In the (p, n) direction the strength distribution 
is dominated by the GT component up to about 20 MeV, 
although non-negligible contributions from IVSM com- 
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FIG. 12: (Color online) Nuclear matrix element for the 2^13(5 
decay of as a function of the coupling strengths Xph 

(a) and Kp^ (b) for prolate shapes in both ^®°Nd and ^®°Sm. 
The shaded area indicates the experimental range extracted 
from the measured half-life using bare qa = 1.273 (lower line) 
and quenched qa = 3 (upper line) factors. 
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FIG. 13: (Color online) Running sums of the 2i'/3P NME in 
'^^^Cd as a function of the intermediate excitation energy in 
'^^^In for various calculations using different quadrupole de¬ 
formations for parent and daughter nuclei. Thin lines cor¬ 
respond to calculations with the energy denominator Di 0 
while thick lines correspond to D 2 ®. 


ponents are found between 10 and 20 MeV, overlapping 
with the GT strength of the GT giant resonance. Above 
this energy, there is no significant amount of GT strength 
in the calculations. The interference effects are construc¬ 
tive below 20 MeV and destructive above this energy, 
but they are not very important. In the (n,p) direction 
the GT strength is expected to be strongly Pauli blocked 
in nuclei with more neutrons than protons and there¬ 
fore, the measured strength is mostly due to the IVSM 



FIG. 14: (Color online) Same as in Fig. [131 but for the double- 
P decay of ^^^Nd. 
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FIG. 15: (Color online) Running sums of the 2vPP NME in 
^^®Cd as a function of the intermediate excitation energy in 
^^®In for the decay between prolate shapes. The results are 
obtained with the energy denominator D 2 for three different 
Skyrme forces. 


resonance over the whole excitation-energy range. Nev¬ 
ertheless, the strength found in low-lying isolated peaks 
is associated with GT transitions because the continu¬ 
ous tail of the IVSM resonance is very small at these 
energies and is not expected to exhibit any peak. In 
summary, the measured strength in the (p, n) direction 
can be safely assigned to be GT in the low energy range 
below 10 MeV and with some caution between 10 and 
20 MeV. Beyond 20 MeV the strength would be prac¬ 
tically due to IVSM. On the other hand the measured 
strength in the (n,p) direction would be due to IVSM 




































[MeV] 


FIG. 16: (Color online) Same as in Fig. 1151 but for the decay 
of i®°Nd. 


transitions, except in the low-lying excitation energy be¬ 
low several MeV, where the isolated peaks observed can 
be attributed to GT strength. This is the reason why we 
plot experimental data in Fig. IHonly up to 5 MeV. 

Calculations for the GT strength distributions in both 
directions B(GT=*=) for the A = 116 isobars are also avail¬ 
able from spherical QRPA calculations within a Woods- 
Saxon basis and two-body Bonn-A realistic interactions 
[s^. The results from this reference are quite similar to 
those in Figs. [3] and IH although a direct comparison 
is not easy as they are obtained for particular values of 
both the strength of the pp residual interaction and the 
effective gA value. The B(GT“) strength distribution 
calculated in Ref. is characterized by some marginal 
strength up to an excitation energy of 5 MeV, then a 
small bump develops between 7 and 10 MeV, and fi¬ 
nally the GT resonance appears between 11 and 15 MeV. 
Both the energy location of the bumps and the strength 
contained agree with the features observed in Fig. |31 
Similarly, the B(GT+) distributions in Ref. [s^ and in 
the corresponding Fig. |3]in this work, contain a similar 
amount of strength concentrated in the very low-lying 
excitation energy. 

Figs. O and m show the results for ^®°Nd and ^®°Sm, 
respectively. We compare the QRPA results from SLy4 
for prolate and oblate shapes with the data from (^He,t) 
on ^®°Nd and from (t,^He) on ^®°Sm in Ref. Q. In gen¬ 
eral terms, we reproduce fairly well the global properties 
of the GT strength distributions, including the location 
of the GT resonance and the total strength measured, 
although the experimental GT resonances appear more 
fragmented than the calculated ones. Similar comments 
to the previous case A = 116 relative to the contributions 
from GT and IVSM components apply now to the case 
A = 150. Microscopic QRPA calculations with residual 
interactions obtained from the CD-Bonn potential were 


performed in Ref. Q to estimate the GT and IVSM com¬ 
ponents of the strength B(GT-I-IVSM) measured. From 
this analysis it was also concluded that the strength 
B(GT“-|-IVSM“) below 20 MeV correspond to a large 
extent to GT strength, whereas the measured strength 
B(GT+-|-IVSM+) below 3 MeV could be assigned to the 
GT component. 

In Fig. [7]we show the results for the GT^ strength dis¬ 
tributions obtained with three different Skyrme forces to 
see the sensitivity of these distributions to the nucleon- 
nucleon effective interaction used in the calculations. We 
can see that the profiles obtained with the three forces are 
quite similar. The GT“ resonances appear at close en¬ 
ergies and contain about the same strength for the three 
interactions. The GT’*' strength is very small compared 
to the GT“ and is concentrated for the three forces in 
the excitation energy region below 5 MeV, with a bump 
in ^®°Sm centered at about 15 MeV. 

Figures [51 [51 [TUI and m show in more detail the GT 
strength distributions with SLy4 in the low excitation- 
energy range. In the case of ^^®Gd (Fig. [5]) we have in¬ 
cluded, in addition, the experimental GT strength mea¬ 
sured by ground-state to ground-state electron capture 
on ^^®In [^. This value agrees well with the strength 
obtained with the prolate shape and to a less extent with 
the oblate one, but it is too low compared to the spheri¬ 
cal calculation. Similarly, in the case of ^^®Sn (Fig. [U]) we 
show with an open circle the experimental GT strength 
measured by the ground-state to ground-state /3“ decay 
(^^®In —)> ^^®Sn). One should notice that the B(GT) 
strength from Ref. Q is calibrated to this value. Al¬ 
though a detailed spectroscopy is beyond the capabil¬ 
ities of our model and the isolated transitions are not 
well reproduced by our calculations, the overall agree¬ 
ment with the total strength contained in this reduced 
energy interval, as well as with the profiles of the accu¬ 
mulated strength distributions is reasonable, especially 
for the prolate shapes in both A = 116 partners in Figs. 
[8]and[9l as well as with the prolate profiles in the A = 150 
partners displayed in Figs. [TU] and m 


B. Double-beta decay 

It is well known that the 2 j^/3/3 NMEs are very sensi¬ 
tive to the residual interactions, as well as to differences 
in deformation between initial and final nuclei [isl . 
To illustrate the dependence of the NMEs to the resid¬ 
ual forces, we show in Fig. [TU] the SLy4-QRPA results 
with prolate shapes for the ^®°Nd ^ ^®°Sm for various 
strengths of the coupling constants Xph The 

shaded region corresponds to the experimental NMEs ex¬ 
tracted from the measured 2uP(3 half-lives, using bare 
{gA = 1.273) or quenched {gA = 1) values. We can see 
that the experimental NMEs contained in the shaded re¬ 
gion are reproduced in windows of the parameters that 
include the values used in this work {Xph = 0-09 MeV 
and Kpp = 0.055 MeV in ^®°Nd). 
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In the next figures we show the running sums for the 
2vfifi NMEs. These are the partial contributions to the 
NMEs of all the I’*' states in the intermediate nucleus up 
to a given energy. We see in Fig. [T3]the running sums for 
the 2i//3/3 decay in ^^®Cd for three combinations of the 
initial and final nuclear shapes. The NMEs of crossing 
deformations produce much lower values than the exper¬ 
imental window and are not shown. Similarly, we show 
in Fig. [14] the running sums for the ^®°Nd ^®°Sm 

decay using prolate shapes for both nuclei. Results ob¬ 
tained with the energy denominator Di © are displayed 
with thin lines, whereas results obtained with D 2 (jS]) are 
shown with thick lines. We can see that D 2 denomina¬ 
tors produce larger NMEs than Di. The main difference 
is originated at low excitation energies, where the rel¬ 
ative effect of using shifted energies is enhanced. The 
effect at larger energies is negligible and we get a con¬ 
stant difference between Di and D 2 , which is the differ¬ 
ence accumulated in the first few MeVs. We also observe 
negative contributions coming from the region of the GT 
resonances at about 15 MeV. This feature was studied 
in Ref. and was related to the strength of the pp 
residual interaction. The small contribution to the 2i//3/3 
NMEs from the region of the GT resonance is due to the 
joint effects of large energy denominators in Eq. and 
the mismatch between the branches of the GT“ and GT+ 
excitations in the energy region of the GT resonance. 

In Figs. US] and HSI we have the 2vl3j5 NMEs calcu¬ 
lated with the denominator D 2 for three Skyrme forces 
in ^^®Cd and ^®°Nd, respectively. The results correspond 
to the prolate shapes in both initial and final nuclei. The 
sensitivity of the NMEs to the Skyrme interaction is man¬ 
ifest in the figures, but it is not as significant as the sensi¬ 
tivity found to other aspects such as the residual interac¬ 
tions, the deformation, or even the energy denominators. 

Another issue worth to comment is the treatment of 
the quenching factor of the axial-vector coupling con¬ 
stant Pa- The physical reasons for this quenching have 
been studied elsewhere and are related to the 

role of non-nucleonic degrees of freedom, absent in the 
usual theoretical models, and to the limitations of model 
space, many-nucleon configurations, and deep correla¬ 
tions missing in these calculations. The implications of 
this quenching on the description of single-/3 and double- 
/3 decay observables have been considered recently in sev¬ 
eral works [12, ilbl - fs^l . In those works both the effec¬ 
tive value of PA and the coupling strength of the resid¬ 
ual interaction in the particle-particle channel are con¬ 
sidered free parameters of the calculation. The strik¬ 
ing result found is that very strong quenching values are 
needed to reproduce simultaneously the observations cor¬ 
responding to the 2z//3/3 half-lives and to the single-/? de¬ 
cay branches, namely log ft{EC) values for the ground 
state to ground state decay of the intermediate nucleus 
into the initial one, and the log values for the 

ground state to ground state decay of the intermediate 
nucleus into the final one. Different procedures for ad¬ 
justing those parameters have been explained in Refs. 


32, . The most recent analysis performed in Ref. 

32j| within a QRPA formalism comes to the conclusion 
that considerable quenching is needed, with an average 
value {pa) ~ 0.6 ± 0.2 (g = 0.47 ± 0.16). The need 
of a strong quenching is also required in other theo¬ 
retical approaches, such as different QRPA calculations 
(g = 0.6 — 0.7) [13, IBM-2 models [111, where the 
quenching obtained is pA = 0.41 (g = 0.32) for ^^®Gd 
and PA = 0.35 (g = 0.28) for ^®°Nd, or shell-model cal¬ 
culations that found quenching factors in the interval 
g = 0.45 — 0.74 [Hi for double-/? emitters from “^^Ga up 
to i5°Nd. 

One should note that the QRPA calculations leading 
to the strong quenching that fits the 2vPP NMEs have 
been performed within a spherical formalism neglecting 
possible effects from deformation degrees of freedom. Be¬ 
cause the main effect of deformation is a reduction of the 
NMEs, deformed QRPA calculations shall demand less 
quenching to fit the experiment. Although the goal of 
this paper is not to extract consequences on the quench¬ 
ing needed to fit the experiment, it will be very interest¬ 
ing to explore in the future the implications of deforma¬ 
tion. This study will be part of a more ambitious project 
aiming to extract an effective pA value that fits the decay 
observables taking into account the uncertainties related 
to the theoretical description and including all possible 
2^/?,5-decay candidates. 


IV. SUMMARY AND CONCLUSIONS 

In summary, using a theoretical approach based on 
a deformed HF-l-BCS-l-QRPA calculations with effec¬ 
tive Skyrme interactions, pairing correlations, and spin- 
isospin residual separable forces in the ph and pp chan¬ 
nels, we have studied simultaneously the GT strength dis¬ 
tributions of the double-/? decay partners (^^®Gd, ^^®Sn) 
and (^®°Nd, ^®°Sm), reaching the intermediate nuclei 
^^®In and ^®°Pm, respectively, as well as their 2^/?/? 
NMEs. Our results for the energy distributions of the 
GT strength have been compared to recent data from 
charge-exchange reactions, whereas the calculated 2^/?/? 
NMEs have been compared with the experimental val¬ 
ues extracted from the measured half-lives for these pro¬ 
cesses. 

We have discussed the sensitivity of our results to the 
various ingredients in the theoretical formalism. Namely, 
to the Skyrme force, to the residual interactions, to de¬ 
formation, to effective pA values, and to the treatment 
of the energy denominators of the NMEs. We found dif¬ 
ferent sensitivities to them that have been analyzed and 
discussed. All in all, the method used in this work has 
demonstrated to be well suited to account for the rich va¬ 
riety of experimental information available on the nuclear 
GT response. The global properties of the energy distri¬ 
butions of the GT strength and the 2^/?/? NMEs are rea¬ 
sonably well reproduced, the exception is the detailed de¬ 
scription of the low-lying GT strength distributions that 
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could clearly be improved. It will be interesting in the 
future to extend these calculations by including all the 
double-/? decay candidates and to explore systematically 
the potential of this method. 
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